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Abstract 

A full three-band model for the Cu02 plane with inclusion of all essential inter- 
actions - Cu-0 and 0-0 hopping, repulsion at the copper and oxygen and between 
them - is considered. A general procedure of the low-energy reduction of the primary 
Hamiltonian to the Hamiltonian of the generalized t-t'-J model is developed. An 
important role of the direct 0-0 hopping is discussed. Parameters of the effective 
low-energy model (the hopping integral, the band position and the superexchange 
constant J are calculated. An analysis of the obtained data shows that the experi- 
mental value of J fixes the charge transfer energy A = (cp — e^) in a narrow region 
of energies. 
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I. Introduction 



Since the discovery of high-T^ superconductors, major theoretical efforts were devoted to 
finding the simplest model which would contain all details relevant to superconductivity. 
Anderson first suggested |l| that the two-dimensional single-band t-J model fits well for 
this role. 

Now there is a general agreement that the Cu02 planes are common to all high- 
temperature superconductors. A realistic model for these planes was proposed from the 
first principles by Emery Varma, Schmitt - Rink and Abrahams P], and Gaididei and 
Loktev |]. 

During last years there was a polemic about the equivalency of the t-J model and 
the Emery model [Q] in the low-energy limit. The principal step was made in the work 
of Zhang and Rice where they proposed the idea of the local singlet and idea of the 
Wannier representation for 0-states. The problem of the low-energy reduction has been 
intensively investigated The singlet - triplet effective Hamiltonian was obtained in 

works of Yushankhai and Lovtsov |P] and Shen and Ting [^. In our recent work the 



quantitative comparison of the exact solution for the three-band model with the solution 
for the generalized t-J model was performed. However, the situation seems not completely 
clarified since (1) most of the above-named works were dealing with the unrealistic region 
of parameters t <^ f/^. A, where t is the Cu-0 hopping integral, Ud is the Coulomb 
repulsion at the Cu site, A = (e^ — e^) is the charge-transfer energy, and (2) all above 
named works did not consistently take into account the direct 0-0 hopping. 

In this work we develop a general approach to the consistent low-energy reduction of 
the three-band model. For this we following by Zhang and Rice transform the primary 
Hamiltonian of the model from the terms of the usual oxygen states to the terms of 
symmetrical and antisymmetrical orthonormalized oxygen states on O - clusters. Further 
we introduce a new basis of the local states with the certain number of particles. Since 
only the systems with filling close to unity are of interest , we restrict ourselves to one- and 
two- hole local states. After that, the diagonalization of the local part of the Hamiltonian 
presents no special problems. It should be noted that we consider the direct 0-0 hopping. 
We have found that it does not change crucially the picture of the local states but play 
an important role for the effective hopping parameters. Next, we hold only the lowest 
two-hole local singlet state and consider its transitions over the background of the lowest 
one-hole (spin) states. Thus, we get the t-t'-J model. All other states we take into 
account perturbatively, by applying the canonical transformation. In this way we obtain 
the general form of the second-order corrections to the t-t'-J model and the expression 
for the super exchange constant J . 

At the beginning, we investigate three limiting cases of the complete three-band model. 
Because the Gilbert spaces in these limiting cases are restricted, mathematical treatment 
is simplified. Further, we investigate the general case of the three-band model in the 
region of parameters where the charge-transfer insulator |n] is the ground state of the 
undoped system. We conclude that the t-J model is vahd for the doped charge-transfer 
insulator, i.e. corrections to it are small. Corrections to the hopping integrals at the 
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second and other neighbors are not relatively small, and thus the simple t-t'-J model does 
not follow from the three-band model. Since the experimental value of the superexchange 
constant J is known very well, we determine the value of the charge-transfer energy A 
through J. We have found that a narrow region of energies for A is suitable. 

We also consider the repulsion at the O site Up and the repulsion between the Cu 
and O ions Vp^ terms in the framework of our approach. Their influence on the values of 
A and the first hopping integral ti is discussed. Thus, we perform a consistent and full 
consideration of the three-band model. 

In Sec. II we represent the three-band Hamiltonian in terms of new oxygen states. 
In Sec. Ill we diagonalize the local Hamiltonian and rewrite it in terms of the Hubbard 
operators. In Sec. IV the low-energy reduction procedure is developed. In Sec. V three 
limiting cases of the three-band model are considered. In Sec. VI the quantitative analysis 
is performed. In Sec. VII the Up and Ypd terms are considered. Section VIII presents our 
conclusions. In Appendices A and B the explicit form of the required matrix elements are 
presented. In Appendix C the three-hole states are considered. In Appendix D the shifts 
of the energies of local states connected with Up and Vpd terms are presented. 



II. Three-band Hamiltonian 

The three-band model studied in this paper was originally proposed by Emery 0. It 
is given by the following Hamiltonian 

H =e,Y.< + ^pll + U,Y. « + AH + H' (2.1) 

l,a m,a I 

AH = f/pE<T<i + ^prf E (2-2) 

m {lm),a,l3 

where / and m denote summation over the Cu- and 0-sites respectively, nf^ = dl^dia, 
^ma — vtaPia^ {dio) creatcs (annihilates) the Cu (3d2,2_j^2) hole, pf^ ipia) creates (anni- 
hilates) an O (jixi Py) hole. U^ and Up are intrasite Coulomb repulsion at the copper and 
oxygen respectively, Ypd is the intersite repulsion between the nearest Cu and O holes. 
The hybridization Hamiltonian H' includes Cu-0 and 0-0 hopping terms: 

H' =tY. {dlVma + H.C.) 

{lm)a 

+ H.C.), (2.3) 

{mm') a 



where {l,m) denotes the nearest-neighbor Cu and O sites and {mm') denotes the nearest- 
neighbor O sites. The quantities t and tp are positive. In Eq. ( |2.3| ) the sign convention 
corresponds to the change of the signs of the operators at the odd sites, which corre- 
sponds to the shift of quasimomentum space by (vr/a, n/a). As was proposed in |T^, such 
transformation makes all 0-0 hopping constants of the same sign and negative. 
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The three-band model (|2.1|) , (|2.3| ) and its modifications have been studied recently by 
many authors |T9|, |20|. One of the simplifications used in all analytical works was 
neglecting the repulsion at the oxygen and between the nearest - neighbor copper and 
oxygen. To take into account the main effect, we will firstly consider the model - 
( p.3|) without Up and Vpd terms and later include these terms as the perturbation. 

The Fourier - transformed Hamiltonian (|2.1|) - (|2.3 ) with the above named restriction 
is given by 

+Ud E ^T^k.T^ai^^i , (2.4) 

ki+k2=k3+k4 



and 



H' = 2t'^{d^^[pkaxCos{k^/2) + p^ay cos{ky/2)] + H.c.} 



ko 



-tp E{2[cos((/c^ + ky)/2) + cos((/c^ - ky)/2)]p^^.^p^ay + H.c.} (2.5) 



ka 



where 



dk = J2diexp(-ikvi),pi^^^y)= J2 Pmexp(-ikrm) , (2.6) 

I mGx{y) 

m G {xiy)} denotes the x[y) 0-sublattices. Hereafter, the lattice constant is equal to 
unity. It is reasonable to introduce the symmetrical and antisymmetrical orthonormalized 
p-operators combination as 

gk = (cos(/c^/2)pkx + cos(A;j^/2)pky)(l + 7k)~^^^ 

«7k = (-cos(/cJ//2)pkx + cos(A;^./2)pk^/)(l + 7k)~^^^ 

where 7k=(l/2)(cos(k2:a)+ cos(^a)). 

Thus, the pd- hybridization term in H' includes only the symmetrical state In 
terms of these operators ( |2.7] ) we rewrite the Hamiltonian (p.4|), ( p.5| ) as: 

Ho = E "^ka^ka + JK^ka^ka + Q'kaQ'ka) 
k,o; k,a 

+Ud E c?kiT^k2Tf^k3i'^k4i , (2.7) 

ki+k2=k3+k4 

H' =2tEM4a?ka + i^.C.} 
ka 

-tp E{/^k[gka?ka - gka^ka] + MlL^^c + ^iL^ka]} 
ka 

were quantities Ak, /ik, J^k are 

Ak =(l+7k)'/^ 

/Xk =8cos2(A;,./2)cos2(A;,/2)(l + 7k)-^/' , 
z/k = 4 cos(A;^/2) cos(A;j^/2)(cos^(A;a;/2) 



-cos2(V2))(l + 7k)-'/' 
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It is not convenient to treat this Hamiltonian in the k-space due to the strong hole 
repulsion at the copper sites or large value of the parameter in this Hamiltonian. 
For taking into account this term precisely, we return to lattice representation for our 
Hamiltonian. After inverse Fourier transformation, Hamiltonian ( |2.71 ) has the form: 

Ho =edJ2 dladla + tp J2i(lla(lla + ^laQlc.) 
l,a I, a 

H' =2tJ2 Xu'{dlq,^ + H.c.) (2.8) 

{ll')a 

-tp J2 {f^ll'illaQl'a - QlaQl'a) 
{U')a 

+ ^w{(ltA'a + H.C)] , 

where 

{(IhQi) = Z!(^k,gk)exp(-ik(ri)) , 

k 

{A,/x,z/V = {A,/i,z/}(l-l') (2.9) 
= Hi-^' exp(-iik(l - 1')) , 

k 

where summation over k is produced over the Brillouin zone, coefficients A,/i, v decrease 
rapidly with increasing |1 — The values of A, /i and v for small |1 — 1'| are given in Table 
I. It is easy to get that z/q = z/qo = and Vn,n = z^(^x + ray) = 0. 

One can divide Hamiltonian (|2.8| ) into the local and hopping parts 



Hloc - X! diadla + (Cp - l^otp) ^ QlaQlo 
I, a l,a 



La 



+UdJ2d+di^dl-^dii (2.10) 
+2tXoY.{dlqia + H.c.) , 

ll'a 

Hhop = 2t XI ^ll'i^taQl'a + H.C.) 
ll'a 

-tpJ2{f^ll'i(lta(ll'a - qta^'c) 
ll'a 

+ l'll'{.(lL^l'c. + H.C.)} , 

hereafter sum over /, /' means that / 7^ Now, one can discuss the reasons for transfor- 
mation from the Hamiltonian (|2.1| ) - (|2.3| ) to ( p.lOj ) . The copper rf-state hybridizes only 
with the symmetrical oxygen g-state 0. Taking into account the direct 0-0 hopping 
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does not change this picture. After this separation out of the strongly interacting d and 
q local states, one can find the full energy spectrum of the one- or two-hole local states 
and solve the problem of the low-energy two-hole state at the background of the one-hole 
states (spins) consistently. The direct 0-0 hopping only slightly shifts the energies of the 
oxygen states with the opposite symmetry. But its contribution to the effective hopping 
parameters will be very important (see Sec. V). The transformed Hamiltonian ( 2.101) with 
the definitions (|2.9| ) is equivalent to the Hamiltonian (|2.1|) , ( p.3|) without the additional 
Coulomb terms Up and Vpd and describes the three-band model exactly. 



III. Diagonalization of the local Hamiltonian 

In this paper we study the model relating to high-temperature superconductors, which 
are the systems with a nearly half-filled band. Thus, we shall concentrate our attention on 
the case of one hole over unit filling in the framework of the three-band model. Namely, 
we determine the space of one- and two-hole local states of the Cu02 plane. A primary 
set of states in the one-hole sector is 

K)=dM, |g,) = g+|0), 

\qa)^qM- (3.1) 
The two-hole singlet sector has a set of states 

\^) ^ 4410) , \v) ^ g+g+|0) , \x) = {1/V2){d^qt - d^q^M , 

1^) ^ g+g+|0) ' \X) = (l/v^)(4^t - 49t)|0) ' (3-2) 
\fl) = {l/V2Kq+q+ -q+q+)\0) , 

The two-hole triplet sector is 

|r) ^ (4g+,4g+,(l/V2)(4g+ + 4g+)) |0) , 

|f) = (4^4, 4g+, (l/^/2)(44 + 4g4)) |0) , (3.3) 



where |0) is the vacuum state (state without holes on a cluster). Hence, we have three 
one-hole and nine two-hole states. As will be shown below, some of these states do not 
play any role in a low-energy model. In this part of the work we shall consider only the 
local part Hioc of the Hamiltonian ( p.lO| ). It is convenient to express this Hamiltonian in 



terms of the Hubbard operators rather than in terms of the usual Fermi operators. Such 
representation of the Hamiltonian ( |2.10| ) in the one- and two-hole sector in terms of states 
(13-pD has the form: 

HL = E{e^^r" + i^P - /iotp)Xf"''" 

La 
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^loc 



+ (e, + /iotp)Xr" 
+2tAo(Xf"^" + i7.c.)} , 
= Y: { [iUd + 2e,)Xf^ + 2(6, - fiot,)X'i 



+2(ep + i^otp)Xr + \{ed + + l^otp)Xp + e^X^ 



(3.4) 



+2t\o{X^'' + H.c.)\+ie, + ep-fiotp) 

m=±l,0 



m 
I 

Tmrm 



m=±l,0 



+ epXf'"^'" + 2tAo(Xf'"™ + H.c.) 



where the upper index 1,2 marks one- and two-hole sector of the local Hamiltonian, Xf'' 
is the Hubbard operator at the site /: 



X, 



ab 



\a l){bl\ , a, b = {da, qa, qa, V^; Xy Xi Vi ^"^i C"^) fm}, 



(3.5) 



where index a = ±1/2 in is a spin projection, m = ±1,0 in denotes triplet 
components. 

The representation of Hioc ( p.lO ) in terms of the Hubbard operators Eq. (|3.4|) allows 
us to solve simply the one-site problem. The diagonalization of Hl^^, Hf^^ is performed at 
each site independently. After the diagonalization, Hl^^ is given by 



l,a 

±(ep±/iotp)Xr"} , 



(3.6) 



where 



e/,g = -Ai tRi , Ai = (A ± i2otp)/2 , 

R, = (A^ + AfXlY^^ A = (A - fiotp)/2, A = - 

\fa) and \ga) are the lower and higher one-hole cluster states 

I/a) = U\da) - V\qa) , 
\ga) = V\da) + U\qa) , 
U={{Ri + A)/2R,y/\ 
V = {{R, - A)/2R,y/' , 



(3.7) 



(3.8) 



Our approach gives us the reasons to assume that in the low-energy limit at unit filling 
the background of the Cu02 plane consists of the lowest |/Q;Z)-states (|3.8|) at each cluster. 
Such state is a linear combination of the ci-hole state and the orthonormalized symmetrical 
g-hole state at the nearest oxygens. Virtual transitions of holes in the |/a)-states (spins) 
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at the neighbor sites and back provide the antiferromagnetic type of interaction between 
these spins. 

In Hf^^ (|3.4| ) there are different sectors, separated by the square brackets, which are 
diagonahzed independently. The first sector includes the three singlets dd, qq, and dq. One 
can prove that diagonalization of this sector gives the lowest singlet with the eigenenergy 
well below than the eigenenergies of other states. Thus, the diagonahzed two-hole part of 
the Hamiltonian (|3.4|) is 

HL = E \^E_Xr + E E.Xr^ , (3.9) 

where E_ is the energy of the lowest singlet |c), Ey are the energies of the upper states, 
y is the set of the two- hole local states without |c). For the analytical expression of the 
eigenenergies and for the set of eigenstates see Sections V, VI. After diagonalization of 
the local part ( p.4|) of the full Hamiltonian ( 2.101 ), one can introduce the nondiagonal 



Hubbard operators and rewrite the hopping Hamiltonian ( p.lO| ) in terms of transitions 
between local eigenstates at different sites. 



IV. Low-energy reduction 

A. Zero order 

As was noted, we assume that at unit filling there are holes in the lowest |/a;)-states 
(spins) at each cluster. We will demonstrate below that these spins interact antiferro- 
magnetically. A general expression for the superexchange constant J is derived in Section 
IV.D. 

For the case of one hole over unit filling we take into account only the lowest singlet 
and its transitions from site to site. All upper two-hole states and transitions to them are 
projected out. Such procedure gives the t-t'-J model with hopping at all neighbors in the 
explicit form: 

H,^, = E.Y. Xr + E tu'Xl^^'Xf" , (4.1) 

/ IVa 

where |c) is the lowest singlet, \fa) is the lowest one-hole state, tw are the effective 
hopping parameters from / to /' site. In more usual terms this Hamiltonian (f4.1| ) may be 
written as 

H,_,, = E_ 5:(1 - + 4) + Y: tu'-clci'a 

I W,a 
Cla = Qa(l - fll-a) , = (Qa)+ , hi = {cj ■ q) , ^ = C'^C^C'J^Cl. (4.2) 

Here q"^, Qq, are the electron creation and annihilation operators. 

In the next part of the work we will consider the second-order perturbation theory, 
which gives us a criterion of validity of the t-t'-J model ([4.1|) . Thus, if corrections to 
the energy and the hopping integrals for the lowest singlet from virtual transitions to the 
upper states are small, one can argue that this model ( |4.1j ) is valid. 
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B. SchriefFer- Wolff transformation 



By diagonalization of the local part of the primary Hamiltonian (|3.4| ) we obtain a set of 
the local states. For the low-energy processes one can consider only the lowest states and 



include all others perturbatively. We use the Schrieffer-Wollf (S-W) transformation |T7 



H = exp {-S)Hexp (S) , S+ = -S. (4.3) 

for getting the second-order correction to the pure t-t'-J model. On this way we use 
the smallness of the {A,/i, z/};;/ constants for |1 — l'| ^ 0. The first-order generator of 
transformation 5* and second-order correction are given by 

[Ho,S] = ~H', 6H' = {1/2)[H',S], (4.4) 



For the case of one hole over unit filling Hq = Hf^^ + H}^^ + Ht^ti (|3.6| , |3.9| . |4l| ) and 



H' includes all terms of H^op ( |2.10| ) which are relevant to the transition between the 
lowest singlet and all of the upper two-holes states. In terms of the nondiagonal Hubbard 
operators the Hamiltonian H' is given by 

Hhop =T.T.' F/f,.' {xf''''Xl^- + H.C.} 

ll'af} y 

+ E Fl?;y{xf-Xf^ + H.c.] , (4.5) 

ll'af3,y 

where y is the set of the two-hole states, the prime in the sum ( [4.5| ) denotes absence 
of contribution of the lowest singlet state; fa and gP are the one- holes states ( ^Tsl) . 
The second term in Eq.( [4.5| ) is essential only for correction to the energy of the c- 
singlet. The first term in Eq. (|4.5| ) provides corrections both to the energy and to the 
hopping integrals. Fj^f^^ are the matrix elements of transitions between \d') ® and 
\xj3l') ® \yl) states (x = f,g). Their explicit form is presented in Appendix A. 
In the case of unit filling, Hq for Eq. ([4.4|) is 



Ho = H'^oc + H^c 

and H' is 

H' = E Dy,^,{xy'f-Xr' + H.c.} , (4.6) 

ll'al3,y 

where |0) is the state without holes, D\ii are the matrix elements of transition of two 
holes at different sites into the two-holes state y : \fal)®\f (31') =^ \Ol)^\yl') (see Appendix 
B). Corresponding generators of the S-W transformation are given in Appendices A,B. 

C. Second-order corrections 

Here we derive a general form of corrections to the t-t'-J model. By applying the S- 
W transformation (O), (|4.4| ) to the Hamiltonian H' (4^) one can get correction to the 
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Hamiltonian of c-singlet 
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with 



5He = -Y.MivX^'Ni, (4.7) 
w 



|p/,y|2 |p9'y|2 



^j/ - F- y Ey - + ^9 ~ (^f' 

and Ni = Xp + X^K 

Corrections to the hopping Hamiltonian have the form [IT 



Inl'aP y 

+ZyXfr^'X';;f^{a^pS^)} , (4.9) 

where 

S« = (l/2)cra/3-^r^ , 

here Xy = 1/2 for singlets and = 3/4 for triplets; Zy = 1 for singlets and Zy = —1/2 for 
triplets 

Tlv = -^^^ . (4.10) 

In Eqs. (|4.8|) , ( |4.10| ) we express parameters M and T through the matrix elements without 
spin structure. They are connected with the primary matrix element as 

Fa^,S = , (4.11) 

for |/=singlet and 

F,^^=(2/?)MF,^'^ (4.12) 

for y= triplet, where [ma/3] = (l/2a|l/2/?, Im) are the Clebsh-Cordon coefficients , m = 
±1,0 are spin-1 projections, a,P = ±1/2 are the spin-1 projections, (3 = —j3. One can 
see that the corrections ( [4.7|) , ([4.9| ) to the t-t'-J Hamiltonian ([4.1| ) depend on the filling 
(A^-term in 5He and 5Ht) and magnetic order (S term in 5Ht ). Such terms can not be 
expressed through the simple direct hopping. 

In Section V we will analyze quantitatively the relative magnitudes of the second-order 
corrections. Validity of the t-J model, as well as correctness of inclusion of hopping at 
the next neighbors are well checked by this analysis. 
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D. Superexchange interaction 



In our approach the superexchange interaction arises in the second order of the per- 
turbation theory over hopping of holes at neighboring cluster and back. The situation 
is similar to calculation of the superexchange interaction in the simple Hubbard model 
(4t^/[/). Of course, in the case of small t our result for the superexchange constant must 
be proportional to as was calculated in earlier works @, |l9l in the fourth order of 
perturbation theory. From Eq. ( [4.6| ) one can get 

5H^ = Hj = J2{JuSiS, + Yu'NiNv) . (4.13) 



Such form ( 4.13 ) is more general than in the t-J model. Firstly, there are interac- 



tions between all pairs of spins. Due to the rapid decrease of the constant {A,/i, i/}, 
{J(ii'}23 ~ 10~^J(ii/)), one can omit all next-nearest neighbor terms in Eq. ( [4.13| ) and get 
the Heisenberg term of the usual t-J Hamiltonian. Secondly, Yu/ ^ — (1/4)7^^/ because 
the hole may virtually hop into triplet states and back, i.e. 



Yi, =-E^. , . (4-14) 

with Xy = 2 and Zy = 1/2 for y =singlet, Xy = —1 and Zy = +3/4 for ?/=triplet. Dfi, are 
connected with the matrix elements D^,^^ as follows 

Al',«/3 = (2«)5.;3A1' (4.15) 

for ?/=singlet, 

= [^o^mv (4.16) 

for y =triplet. One can check that in the limit t ^ the reduced matrix element Z^f^, is 
proportional to t and the superexchange constant Jw has contribution proportional to t^. 
But such contributions from the singlet and the triplet two-hole states are cancelled, and 
we get usual result which is proportional to t^. 



V. Limiting cases 

In the previous parts of the work we consistently reformulated the three-band model 
and passed from the terms of the primary Fermi operators of the hole at the copper and 
at the oxygens to the one- and two-holes states Hubbard operators on the cluster. For the 
problem of one hole over unit filling we established the general form of the low-energy the 
t-t'-J model ( [4.1|) , the general form of the corrections to it ([4.7|) , ( [4.9|) , and the general 



form of the spin interaction at unit filling ( [4.13| ). Now we will consider three limiting 



cases, all with special constraints on the Gilbert space of the problem. Such constraints 
simplify the mathematical treatments and provide the possibility to get results in the 
analytical form. 
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A. The case of A < [/^ 

This case was considered by Lovtsov and Yushankai We include into this problem 
direct 0-0 hopping and make consideration more complete. 

Such limit (U^ =^ oo) pushes up in energy the state with double occupation of the cop- 
per. Thus, one can ignore |'?/')-singlet state (p. 21). After such simplification diagonalization 



of the two-particle sector of the Hamiltonian (|3.4|) presents no problems 



'bb 



+2fiotpXr + I + 1 (5-1) 

m=±l,0 

+E+ I X^^ + I - (A + A^otp) E ^™ ™ { ' 

y m=±l,0 / m=±l,0 

where c, h are the lowest and highest singlets oi q — d sector, c, h are the lowest and highest 
q — d singlets, tl, t2 are the lowest and highest q — d triplets, r and (p are determined in 
Eqs. ( |3.2| ), (p.3|). All energies are counted from ep (hereafter = 0). E± = — A2 ± R, 
A2 = 1/2(A + 3/iotp), = (A2 + 8t2A2)i/2^ E± = -A ± A = (A - /iotp)/2. New 
eigenstates are connected with the primary set of states ( |3.2|) , ( |3.3| ), as 



|c) = -a\ifi) + b\x) , |6) = 6|</^) + a|x) , 

6 = ((i? + A)/2i?)^/2 ^ a = ((^ _ A)/2i?)V2 ^ (5 2) 

\5) = U\x)-V\fi) , |6) = + t/|r/) , 

\tl) = U\f)-V\0 , |t2) = r|f) + [/|C) , (5.3) 



U and are determined in Eq. (|3.8| ) The energies of the local states at U(i ^ 00 ( |5.1| ) for 
values of the hopping parameters t = 1.4 eV, tp = 0.7 eV and A = 3.66 eV are given in 
Table II. 

Next, we follow the scheme of Sec. III. Explicit forms of parameters of the t-t'-J model 

are 

E^ = E_ = -{A + 3fiotp)/2- R 

tw = 2t\u'Vh{hU + V2Va) + tp^iu'{V2aV + Ubf /2 , (5.4) 

Hopping parameters for the first four neighbors at values of parameters mentioned above 
are given in Table II. One can see the important role of the direct 0-0 hopping. Its 
contribution to the full amplitude of hopping at the nearest neighbor 50%. Explicit 
form of the matrix elements Ff^f and Df^, one can get from Eqs. ( p.2| ), ( [4.5| ), ( |4.6| ), and 
( p.lOp . They are given in Appendices A and B. Relative magnitudes of corrections to the 
hopping and energy are presented in Tab. II. Corrections to the energy and to the first 
hopping parameter are 4.1% and 3.1%, respectively. Thus, the t-J model is valid for this 
case. It seems that inclusion of hopping at the second and third neighbor to the t-t'-J 
model is justified because corrections to them are small enough. 
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B. The case (U^ - A) <CUd 
This case was considered by the authors 



lOl without the 0-0 direct hopping. In Ref. 



T0| we assumed that at unit filhng there is one hole (spin) per copper site. In terms of 



the present work we ignored the admixture of the |g)-state to |(i)-state. We also assumed 
that there are no doubly occupied oxygen degrees of freedom. Both assumptions are valid 
over the parameter ~ t/Ud- Such constraints lead to U =1, V=0 (see Eq. (|3.8|) ) and rather 
simple form of the diagonalized Hamiltonian Hf^^ 




+ E+X^' - {5 + fiotp 



E 

i=±l,0 



Tin T'tn 



m=±l,0 

where E± = -Ai ± i?, Ai = l/2{Ud - 3A - /xotp), i? = (A^ + St^Xlf/"^, A 
/Xoip)/2, and 

|c) = -a\iP) + h\x) , \h) = h\iP) + a\x) , 
b={{R + K)/2Rfl^ , a={{R- K)/2Rfl^ . 

The explicit form of parameters of the t-t'-J model is following 

E, =E^ = -{Ud-3A-i2otp)/2-R 
til' = 2tXii'V2ba + tpfiii>b^/2 . 



{Ud-A 



(5.5) 



(5.6) 



The spectrum of the local energies and first four hopping integrals at = 8 eV and 
A = 3.65 eV are given in Table III. One can see that direct 0-0 hopping is less important 
than for the case A. It is evident, because part of the oxygen degrees of freedom is excluded 
from consideration. Expressions for the matrix elements of transitions to the excited states 
Fj^iF , Dfi, one can easily get from (|575|) , ( ^.5|) , and (|2.10|) . They are given in Appendices 
A and B. A correction to the energy is not so small as for the case A. (see Table III). It is 
close to 13%, a correction to the first hopping is close to 4.0%. Thus, the t-J model may 
be valid. Inclusion of hopping at next neighbors requires consideration of corrections in 
the form ( |4.9| ) because they are not relatively small. They have a complicated structure. 
One can see a discrepancy in our approach. We assume that (U^ — A) <^Ud, i.e. A ~Ud, 
while the realistic region for A (for which J ~ 126 meV) is close to 3.6 eV. Such situation 
is explained by the importance of the oxygen states for the three-band model. 



C. The special case of the complete model 

Here we consider the complete three-band model (|3.4|) without any constraints on the 
Gilbert space. Let us put (U^ — 2A) = —2fiotp. Such choice of the parameters leads to 
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accidental degeneracy of the energies of the qq and dd singlets. The q — d singlets sector 
of the Hamiltonian (13.41) is written as 



I 

-(A + /iotp)X,^^ + 2V2t\o{Xf + Xr + H.c.)] + . . . , (5.7) 

where = E^^ = Eq = Ud — 2A = —2fiotp. Simplification in this case is as follows: linear 
combination of \ip) and \(p) ( p.2|) singlets 



|c2) = (l/v^)(|^)-|v,)) (5.^ 



does not hybridize with the qd singlet \x) ( |3.2| ). After diagonalization of iff^^ we have: 
I 

+2iiotpXr + E_ix['+ J2 Xf™*i™) (5.9) 

\ m=±l,0 / 

+E+ Ixf'^ + J2 '^"^ I - (A + i^otp) E ™ 

y m=±l,0 / m=±l,0 

where E± = -Ag ± R, Eq = -2jiotp, As = ^7^/4 + 2/iotp, ^± = -A ± A = f/rf/4, 
i?i = (A^ + At'^Xiy/'^, and i? = (A^ + let^A^)^/^ Thus, the real hybridization parameter 
is 16t\o/Ud ! It means that the perturbation approach over t/Ud 0, is unjustified for 
this system. The states |c), \b), \tl) and \t2) are defined in Eq. (|5.2| ) 



|cl) = -a(|v9) + m/V2 + , |c3) = 6(lv^) + m/V2 + a|x) , 

b = {{R + A)/2Ry/^ , a = {{R- K)/2Rf'^ , (5.10) 

Since in this case A =11^/2 + /iotp, one can determine an adequate value of at fixed t, tp 
and J. In the case being considered, the experimental value of J is achieved at = 8.13 
eV and A = 5.1 eV. Expressions for the parameters of the t-t'-J model are 

E,i =E_ = -{Ud/A + 2/^otp) - R. 

tw = 2t\u' {UV + ha) + tp^in> {aV + hUf /2 . (5.11) 

The local energies and hopping parameters for this special case are given in Table IV. 

The matrix F^f and Df^, elements which are required for the derivation of corrections 
to the energy E_, the hopping constants tui and the superexchange constant J are given 
in Appendices A, B. As in the first limiting case, 0-0 hopping plays an important role 
for the magnitude of effective hopping of the lowest cl-singlet. Correction to the energy 
is close to 4.6%, correction to the first hopping constant is close to 2.5%. Thus, the 
t-J model is valid with the same precision. Absolute magnitude of the second hopping 
integral (hopping at the next-nearest neighbors) is small due to the partial compensation 
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of the amplitudes of Cu-0 and 0-0 hopping (see Eq. and Table IV). As a result, 

correction to them is rather large ~ 78%, and the t-t'-J model is not valid. Relative 
magnitude of terms with transition at the next neighbors (farther than nearest) is of the 
order of 10%. This is also the parameter of accuracy of the t-J model as the low-energy 
limit of the three-band model. 



VI. The general case 

In this part of the work we consider reduction of the three-band model to the effective 
t-t'-J model in the general region of parameters when the charge transfer insulator is the 
ground state of undoped system. The general case differs from the special case considered 
above only by the absence of degeneracy of the local \dd) and \qq) two- hole states. The 
local Hamiltonian is defined by Eq.(p.7|) with =11^ — 2 A and = — 2/itp. We stress 
that E^ 7^ E^ in this case. The expression for the diagonalized Hamiltonian Hf^^ coincides 
completely with Eq.( |5.9| ) but the eigenstates \ci) have the more general form 

\cz) = U,\ij) + V,\y,) + W,\x), (6.1) 

where coefficients Ui,Vi,Wi are determined from the solution of a system of three lin- 
ear equations. The energies E^, Eq, E+ are roots of the corresponding cubic equation. 
Effective hopping parameters are 

tw = 2tXiv{V2UiU - WiV){V2ViV - W^U) + tp^iw{V2ViV - W^Uf /2 , (6.2) 



Other matrix elements which are required for the derivation of the constants F^l^ ( [4.5| 
(IJ), ( CT) and Dl, (^), (^^), (^TeJ) are presented in Appendices A,B. 



Parameters of the effective t-t'-J model are calculated at the following values of pa- 
rameters of the three-band model: t = 1.4 eV, tp = 0.7 eV, = 8 eV and A = 5.1 eV 
according to different band calculations |jl3|, |14|, ^ and the detailed analysis of experiment 



In Table V the magnitudes of the hopping parameters at the first four neighbors are 
given. The important role of 0-0 hopping is shown: for the first hopping parameter it is 
close to 30%, for others is close to 50%. Relative magnitudes of corrections to the energy 
and first hopping are 3% and 2.5% respectively. Thus, the t-J model is valid with the 
same precision. Due to compensation of contribution of the p — d and p — p hopping to 
the effective hopping integral at the second neighbors (they have the opposite sign), the 
correction to it is not small (82%). Thus, the model with transitions at the neighbors 
farther than nearest must include the terms of the form ([4.9|) , which depend on filling and 
the spin state of neighbor sites. 

Applicability of a perturbation scheme in the realistic region of parameters of the 
three-band model is provided by smallness of the ratio of the effective hopping parameters 
between different local states to the energy gap between them. This ratio is of the order 
of 10% which gives the accuracy of the t-J model for the Cu02 plane. Relative magnitude 
of the t' -terms (hopping at the next neighbors) in the t-t'-J model is also of the order of 
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10%. In Table V the fundamental ratio ti/ J = 4.2 for = 8 eV is presented. It only 
weakly depends on and A. This value is slightly larger than the generally accepted 
t/J^ 3. 

Different parameters of the three-band model are known with different accuracy. No- 
tice that the value of the superexchange constant for La2Cu04 is known with high accu- 
racy which imposes restrictions on the values of the parameters of the three-band model. 
Therefore one can try to determine the value of the least known parameter A as a function 
of other parameters of the model. At fixed values of t and tp presented above and J = 126 
meV we get A as a function of in the physically reasonable region: A < < 12 eV 
in Fig. 1. One can see that the reasonable region of A lies between 4.75 eV and 5.75 
eV which is in complete agreement with the band calculations of Sushkov and Flambaum 



1 15 1 but larger than the generally accepted value A = 2.75 — 3.75 eV |22 



VII. Inclusion of V^^ and terms 

Here we consider the role of the intrasite repulsion at the oxygens Up and the intersite 
repulsion Vpd- We suppose that the main effect is taken into account ( strong hybridiza- 
tion d- and g-states), so that one can take the above named terms as the perturbation. 
Therefore, we do not consider the contribution of these terms to the second-order correc- 
tions to the energies and the hopping integrals of the form ( |4.7] ), ( [4.91 ). Our aim is to 
get the renormalization of the first hopping integral and the charge-transfer energy due 
to the Hamiltonian AH (|2.2| ). Using the representation ( p. 71) for qi and qi operators one 
can easily get 

AHpd =Vpd E < {/^^(li - 1> I2 - l)g/t/3%,/3 

I,ll,l2,a,l3 

+ - 1, 12 - l)[gt/3*2,/3 + H.C.] (7.1) 

Since we are interested in the renormalization of the lowest singlet parameters, we can 
omit all terms with antisymmetrical oxygen operators q. Similar procedure for Up - term 
leads to 

AH = Vpd E fi,iuhnfaltp1h,P 

+Up J2 Kiuh,U(l^(lhA)iltl1i3d) ' (7-2) 

where fi^i^^i^ = /^^(li - 1, 12 - 1), /i«,h,/2,«3 = ^(^i -hh- 1, 13 - 1), and 

fUh - 1, 12 - 1) = E /k,k' exp(^k(li - 1) - zk'ih - 1)) 

k,k' 

/i(li- 1,12-1,13-1) = E ^k,k',k" (7.3) 

k,k',k" 

exp(zk(li - 1) - zk'(l2 - 1) + zk"(l3 - 1)) 
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with 

fk,k' = 2{cos(A;,/2) cos(A;;/2) cos((A;, - k'J/2) + (x ^ y)}(l + 7k)~'/'(l + 7k')"'/' 
hk,k',k" = 2{cos(A;,/2) cos(A;;/2) cos(A;;V2) cos((fc, + K- k'i;)/2) (7.4) 
+ (x ^ y)}(l + 7k)-'/'(l + 7k')"'/'(l + 7k")"'/'(l + 7k+k'-k")"'/' • 

Thus, we have the four-fermion terms with the comphcated relation between many states 
at different sites. All later arguments are based on a very rapid decrease of the constant 
fw and hi^i^^i^^i^ with |ln — 1|. Our calculation gives /o = /(0,0) = 0.9180, fi = /(1,0) = 
/(0,1) = 0.1343, and Sf = Ei.^i.^i fh,h,h = 0-0092; ho = /i(0, 0, 0) = 0.29 and h = 
/i(l,0,0) = /i(0, 1,0) = /i(0,0, 1) = 0.0096. Therefore, one can turn from ([Tj) to the 
effective Hamiltonian 

AHlil =V,,h E + i^.c.} (7.5) 

+Uphi E {n1^{qi^iqi',l + H.c.} + n'^iiq^qi',^ + H.c.}} , 

m 

where in AH^J/ fw = f{V - 1, 1' - 1). 

Let us consider the system at unit filling and calculate the shift of the energies of the 
local states due to AH^J/ . At unit filling there is the state \ fa) = U\da) — V\qa) at each 
site. Hence, the shift of the energy of d - hole at the site / will be 

Ae, = Vp,V^ E fw = VpdV\2 - /o) (7.6) 

for the qi state 

Ae, = VpdU\2 - /o) + UpV\l/2 - ho) (7.7) 

numbers 2 and 1/2 are the sums over all I' for fi'^i and hf^i respectively. The coefficients 
U and V are defined in Eq. (|3.^ ) with 



R, = {A' + 4t'\lf\ A = (e,-e,)/2, (7.^ 



Thus, Eqs.( [7.5D , (|7.6| ), (|7.7|) , and are the system of equations, which can be solved 
numerically. After solving we have new values U and V, and consider one hole over unit 
filling. It is easy to get the shift of the local energies of the additional hole. For some 
details see Appendix D. 

The matrix elements of AH^{^ (|7.4]) produce the addition to the first hopping integral. 

Ati = {-fiVpdUWi + V2hiUpVVi){WiU - VWiV) (7.9) 

The expressions for the matrix elements AD^i, which are essential for the constant J ( [4.14|) , 
are given in Appendix B. 
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Different band calculations |jT3|, ^ give the consistent magnitudes for Vp^ and U- 
Vpd = 1.2 eV, Up = 4 eV. The calculations present no special problems. Fig. 2 and Table 
VI present our results. The lower curve in Fig. 2 demonstrates that for reasonable region 
of (5 — 11 eV), at fixed Vpd, Up, t and tp, A varies in the region 2.5 — 3.5 eV, which 
coincides very well with the band calculations ^l| and analysis of experiments 



Different contributions to the first hopping parameter tl are presented in Table VI for 
U^ = 8.0 eV and A = 3.0 eV. Inclusion of AH tends to decrease ti. The ratio ti/ J = 3.38, 
which is close to the results of the other works O, ITS 



VIII. Stability of the system 

There are few conditions which determine stability of the system under study. One 
can take an analogy from the classical Hubbard system: the energy of two electrons at 
different sites ({1/, 1;/}) is lower than the energy of two electrons at the same site and 
an empty site ({2i,0//}) by the energy of the Hubbard intrasite repulsion Uj^. Thus, the 
ground state of the system at unit filling is the state with one electron per site which is 
a dielectric if Uh is sufficiently larger than band width wzt, where z is the number of 
neighbors and w is the constant of the order of unity. 

Since we have reduced our problem to the system of the Hubbard type but with many 
states at each site, we can formulate the stability conditions in terms of energies of these 
states. The conditions will be 

E'^ + E^-2E^>0 (8.1) 

and 

+ -2E^ >0 , (8.2) 

where the indexes n = 1,2,3 denote the states with n holes at site. From this one can 
easily get the charge-transfer gap (an analogue of U^) as the difference between (E'^ + E^) 
and 2E^, where E"^ is the energy of the lowest singlet and E^ is the energy of the lowest 
one-hole state 

AE{2, - 1, 1) = ^_ + - 2^/ = 3.2 eV > (8.3) 

for the above mentioned parameters. 

An interesting quantity is also the difference between {{2i,2ii}) and ({3;, 1/'}) states 

AE{3, 1-2,2) = E^ + Ef - 2E_ = 2.7 eV . (8.4) 

where E^ is the lowest three-hole state (see Appendix C). This energy difference AE{3, 1 — 
2, 2) is also analogue of the energy of the Hubbard repulsion U^- 
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IX. Conclusion 



We conclude by summarizing our results. We have studied the low-energy properties 
of the Cu02 plane near unit filling in the framework of the three-band model. We have 
considered the full three-band model in the form that was put forward by Emery 
in the realistic region of the parameters, without any additional assumptions about the 
smallness of some of them. 

Thus, we have taken into account the direct 0-0 hopping, intrasite repulsion at the 
oxygens Up, and intersite repulsion Vp^, which have not been considered earlier. We have 
presented the consistent approach to the mapping of the three-band model to it low-energy 
limit. Following the idea of Zhang and Rice 0, we expressed the primary Hamiltonian in 
terms of symmetrical and antisymmetrical oxygen states. Next, we turned to the terms 
of the local states with certain number of the particles and transition between them. The 
diagonalization of the local part of the Hamiltonian provided the set of eigenstates. The 
lowest of them is the singlet state with the energy well below than others. We have derived 
the low-energy Hamiltonian for the lowest singlet and its transitions [t-t'-J model), and 
have taken into account all upper states by the special type of the unitary transformation. 

Thus, we obtained the effective single-band Hamiltonian which is essentially the t-J 
model one. Transitions to the next neighbors are not simple hopping due to the important 
role of the correction, which has a complicated structure. 

Our approach allows to establish the quantitative boundary of the validity of the t-J 
model as the low-energy limit of the three-band model, to get the corrections to it in an 
explicit form, and to take into account the transitions at the next neighbors. It is evident, 
that one can determine the value of the charge- transfer energy A from the well defined 
value of J at fixed other parameters of the three-band model. We have established that 
A varies in a narrow region of energies. 

The lowest one-hole state for the undoped Cu02 plane is 



\fa) = U\da) -V\qa) a =T, i • (9.1) 
The lowest two-hole singlet state for the doped Cu02 plane is (|6.1|) 

|cl) = U,\d Ui) + Vi\q Ui) + W,{\d ]qi)-\diq ]))/V2 (9.2) 

where \da) is the state of a hole at the copper with projection a ; \qa) is the state of a 
hole at the symmetrical oxygen which represents the Wannier state formed from the four 
oxygen states around the copper ion; \d ^ d I) , \q '\ q i) and {\d '\ q I) — \d I q t))/-\/2 
are the two-hole singlet states. 

The following set of parameters of the three band model ( p.l|) - ( p.3|) seems the most 



acceptable at present: A = ep — = 3 eV, Ud = 8 eV, f/p = 4 eV, Vpd = 1.2 eV, t = 1.4 
eV, tp = 0.7 eV. 

For this set of parameters of the three band model we have the following values of the 



coefficients f/, V, f/i, Vi, Wi (lOP, ( |9.2| ) which determine the probability of location of 
holes at the copper and the oxygen: U = 0.85, V = 0.52, Ui = -0.38, Vi = -0.64, Wi = 
0.67. 
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The fundamental parameters of the t-J model ti and J are : ti — 0.427 eV, J — 0.126 
eV and ti/J = 0.34. 

The charge transfer gap AE — 3.2 eV and the effective Hubbard repulsion of holes 
Ui/ = 2.7 eV. 



20 



Acknowledgments 



This work was supported, in part, by a Soros Foundation Grant awarded one of us 
(V.I.B.) by the American Physical Society. This work was supported partly by the Soros- 
Akademgorodok Foundation, by the Counsel on Superconductivity of Russian Academy 
of Sciences, Grant No. 90214, and by the Program 'Universities of Russia as the centers 
of fundamental researches'. 



Appendix A: Matrix elements for the singlet tran- 
sitions 

Here we present the explicit form for the matrix elements of the transition of the 
additional hole from the lowest to excited states. Thus, 

Ff^;^' = {ly\{l'xf3\\HHor,\\d')\fal), (A.l) 

where H^op determined in eq. (|2.10|) . A generator of the S — W transformation is 

S = - E E ' {Xf '^"Xf -H.c.} 

Wap y ^- 
p9f^,y 

- E r, _ {Xf-^°x,f'--g.c.} . (A.2) 

ll'al3,y ^y + ^9 

As was noted, the spin structure of the matrix elements can be separated out (see Eqs. 
( [4.11| ), ( [4.12D ). So, for the first limiting case (Sec. V.A) the matrix elements are 

F/iI'' = 2t\ii,V{Uab - ib'^ - a^)V/V2) - tp^iw{VhV2 - Ua)(Uh + VaV2)/2 , 

Flf = 2t\wV'^a - tpfiu>U{Ub + VaV2) /V2 , 

F// =tpUu,{Ub + VaV2)/2 , 

Flf = -tpvw {Ub + VaV2)/V2 , 

Fi^ = Flf = Fff' ^ , (A.3) 

F^if = 2t\w{b{U'^ - V'^)/2 + aUVV2) - tp^iwiVb - V2Ua){Ub + VaV2)/2 , 

F^f = tXiv{ab{U'^ - V"^) - UV{1? - a^)y/2) - tpfmiUa - VbV2){Vb - UaV2)/2 , 

Ff/ = -V2tXivV{b - UVaV2) + tpfiwU{Vb - V2Ua)/V2 , 

Fif =-tpiya>{Vb-UaV2)/2 , 

Fif = -tpuiv {Vb - UaV2)/V2 , 
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For the second limiting case (Sec. V.B.) there are matrix elements only with x — f 

Flf' = -s/2t\ii,{h'' - a") + tppiii,ba/2 , 

f// = 2tXwa + tpHwb/V2 , 

F// = tpUu'b/2 , (A.4) 

F// = tpVu'h/V2 . 

The matrix elements for the third case (Sec. V.C.) 

F^f = tXivKV - 1/2) - tpi^u'V{Ub + Va)/2 , 

Flf = tXirib-" - a^) + Vz,((62 _ ^2^y^ _ ^^2 _ v^)ha)/2 , 

F^/ = y/2t\ir{U^ - V^)a + tpixu,U{Ub + Va)/V2 , 

F// ^-tpPu'{Ub + Va)/2 , 

F/f ^-tpPu,{Ub + Va)/V2 , 

Fif = F// = F/f ^ , (A.5) 

Ff/' = tXu'iU^ - V^) - tpi^uiVb - Ua){Ub + Va)/2 , 

Ff/'^ = t\u'{2UVb -a)- tp^iviVb - Ua)V/2 , 

Fff =tp^iv{Ua-Vhfl2 , 

Ff/ = V2tXwV{2UVa - b) + tpinvUiVb - Ua) /\/2 , 

F?if ^-tpiyu'{Vb-Ua)/2 , 

F/f =-tpUn'{Vb-Ua)/V2 , 

II 11 It 

The matrix elements for the general case (Sec. VI.) are 

f//' = -tXu>[{V2UiU - WiV){V2ViV - WiU) + {y/2ViV - WiU){V2ViV - W^U)] 

-tp^iw{V2V,V - W,U){V2ViV - WiU)/2 , 
Ff/' = -tXi,[iV2U,V + WiU){V2ViV - WiU) - {V2V1U + WiV){V2UiU - WiV)] 

+tpiiuiV2ViU + WiV){V2ViV - WiU)/2 , (A.6) 

where i = 1, 2, 3 is the number of the singlet from q — d sector. Other matrix elements 
are 

F// - 2tXw{U^V^ - U,U^) - tpfiwU{V2ViV - W^U)/V2 , 

F// = tpiyi,{V2ViV - WiU)/2 , 

F/f = -tpiyuiV2V,V - WiU)/V2 , 

Flf' = f// = Fif ^ , (A.7) 

fI/ = -V2tXu,{Wi + (C/i + U^)UVV2) + tpiiii,U{V2V,U + W,V)/V2 , 
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Ff = -tpUu,{V2ViU + WiV)/2 , 
Fif = -tpui, {V2ViU + WiV)/V2 , 



Appendix B: Matrix elements for the virtual tran- 
sitions 

For the case of unit filling we need to determine the matrix elements of the virtual 
transition of a hole at neighboring clusters. Thus, we are interested in the following 
quantities 

Df,^^p = {lym\H,,,\\fal')\f(3l) . (B.l) 

The spin structure can be separated out (see Eqs. (|4.15|) , ( |4.16|) ). A generator of the 
S — W transformation has the form 

S = -T. i'r! , {Xr^'^xr - H.c.} . (B.2) 

ll'al3,y + -^0 ^t/ 

For the first limiting case (Sec.V.A.) 





= V2tXw{b + aUVV2) + tp^iwV{Uh 


+ VaV2)/V2 , 




= V2tXw{a - bUVV2) + tpfiivV{Ua 


- VbV2)/v^ , 










= tpUu'V/V2 , 






= tpUii'V , 






= Df, = D'^, ^ . 





(B.3) 



For the second limiting case (Sec.V.B.) 

D'li, = -V2t\a,h , 

D\i, = -V2tXwa , (B.4) 

DJii = —2t\ii' , 

Df, = DJ, ^ . 

For the third limiting case (Sec.V.C.) 

Df^ = V2t\w{b + 2aUV) + tpfiii'V{Ub + Va)/V2 , 
^iv = tpfiii'V"^/ V2 , 

D'f^, = V2tXiii2bUV - a) - tpiJu>V{Ua - Vb)/V2 , 
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DJ, =2t\u'iU^-V^)+tpfiu'UV , 

Dl, = tpUwV/V2 , (B.5) 

dI = D% = D'if, ^ . 

For the general case (Sec. VI.) matrix elements are 

D'li, = -V2tXu,{V2UV{U, + Vi) - W,) - tpfiu'{V2V^ - W,U)V/V2 , 

Dl, =2tXa'{U'~V^)+t,i2a'UV , 

D% = tpUu>V/V2 , (B.6) 

D*/, = tpUu'V , 

— n'^ — n*^ = n 

— — ^iv — u • 

The additions to the matrix elements from Up and Ypd terms, which are essential for 
calculation of J (Sec. VII) 

AD$,^ = VpJ.UVW, - Uph,V,V^ , 

ADl,^ = Vp,hUV . (B.7) 



Appendix C: Three-hole states 

There are twenty of the three-hole cluster states 















IK) 




1 ai 








and 










^dtq^q:!;\o) . 







(C.l) 



THe hybridization term ( 2.101 ) provides the transitions only between 

1^^) - l<^>') , (C.2) 

and 

|v^i),|$i)^(2a)|A.„..) . 
Four other states {A^^a) are noninteracting and degenerating in energy 

El = -A (C.3) 

counted from ep = 0. 
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By diagonalization of the Hamiltonian (|2.10 ) one can easily get that the eigenfunctions 
1$^) , and \Aa--aa) in the three-holes sector are connected with the q — d sector of 
the two- hole states (|5.9|) . Three- hole eigenfunctions are 



(C.4) 

and their energies 

E^{Ga,^) = E^-^^ = /iotp + , (C.5) 

where Ef = E±^q are the energies of the local singlets ( |6.1|) . Diagonalization of the other 
sectors ( |C.2| ) provides the eigenenergies 

E| = -(A + 3/iog/2± (((A-/iotp)/2)2 + 4t2A2)'^' , (C.6) 

for the (4>^, <l^) sector 

Ei = ~{U, - 3A - 3^otp)/2 ± (((f/d - A + ^lotp)/2f + W'xtf'^ , (C.7) 



and for the (\E'J?,, sector. 

Since we are interested in the lowest three-hole state, we analyze the spectrum of the 
energies ( |(J.5| ) , ( |U.6| ) . For a wide region of and A 

E^-^^ < El , (C.8) 

and all other states lie much higher. Renormalization of these energies due to the Up and 
\pd terms remains all our consideration to be correct. 

For the above mentioned parameters Ud = 8 eV, A = 3 eV, Up = 4 eV, and Ypd = 1.2 
eV one can get 



^1^+9" = -0.5 eV , 
El = -3.73 eV . (C.9) 

Thus, the Hubbard repulsion is 

Uh = E^ + E^ - 2^2 = -3.73 - 4.27 + 2 ■ 5.36 ~ 2.7 eV . (C.IO) 



Appendix D: Shifts of energies of the two- hole states 

Shifts of energies of the primary set of the two-hole states, due to Up and Vpd terms are 

AE^ = 2VpdJ\2 - /o) + Up{ho + V\l/2 - ho)) , 

AE^ =2VpdV\2-h), 

AE^ = 2Vpd + UpV^il/2 - ho)) , (D.l) 

AEr =2Vpd + UpV\l/2-ho)) . 
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TABLE I. 

Coefficients A(l — 1'), //(I — 1'), and vil — 1') as functions of (1 — 1') = rax + my 



n,ni ^n,m ^m,n f^n,m f^m,n 



: 


,0 


0.9581 


1.4567 


0.0 


1 , 


,0 


0.1401 


0.5497 


0.2678 


1 , 


, 1 


-0.0235 


0.2483 


0.0 


2 


. 


-0.0137 


-0.1245 


0.0812 


2 , 


, 1 


0.0069 


-0.0322 


0.0609 


2 , 


, 2 


0.0035 


0.0231 


0.0 



TABLE II. 

Energies of the local states, the first four hopping parameters, the second-order corrections to 
them on a ferromagnetic background, correction to the energy the superexchange constant 
J, the ratio of {h/J) for the case A <C at A = 3.66 eV. 
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one-hole local 


Ef = -5.33 


Eg = -0.65 


Eg = 1.02 




energies, (eV) 










two-hole local 


E_ = -7.38 


E+ = 6.57 


= -4.68 




energies, (eV) 


E- = -4.31 


E+ = 1.67 


E^ = 2.04 




neighbor 


direct hopping 


contribution of 


contribution of 


second-order 


number, 


parameters, 


direct 0-0 hopping. 


Cu-0 hopping, 


corrections, 


n 


tn, (eV) 


(eV) 


tPf, (eV) 


Stn, (%) 


1 


0.43 


0.24 


0.19 


3.1 


2 


0.078 


0.11 


-0.032 


9.6 


3 


-0.074 


-0.055 


-0.019 


UA 


4 


0.005 


-0.014 


0.009 


110. 


correction to 




superexchange 




ratio 


the energy E_, (%) 




constant J, (meV) 




\ti/J\ 


4.1 




126.1 




3.4 
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TABLE III. 

Energies of the local states, the hopping parameters at the first four neighbors, the second- 
order corrections to them on a ferromagnetic background, correction to the energy E_, the 
supcrcxchange constant J, and the ratio ti/J for the case U(i — A ^ at = 8 eV, 
A = 3.65 eV. 



one-hole local — —3.65 

energies, (eV) 



two-hole local 


E_ = -6.63 


E+ = 2.66 


E^ = -4.68 




energies, (eV) 


E^ = -2.63 


Ef = -2.63 






neighbor 


direct hopping 


contribution of 


contribution of 


second-order 


number, 


parameters. 


direct 0-0 hopping, 


Cu-0 hopping, 


corrections. 


n 


in, (eV) 


t^i, (eV) 


t'k', (eV) 


5tn, (%) 


1 


0.38 


0.15 


0.23 


4.6 


2 


0.030 


0.068 


-0.038 


313. 


3 


-0.057 


-0.034 


-0.022 


91. 


4 


0.0022 


-0.0089 


0.011 


165. 


correction to 




superexchange 




ratio 


the energy E_, (%) 




constant J, (meV) 




\ti/J\ 


13.0 




126.6 




3.0 



29 



TABLE IV. 

Energies of the local states, the hopping parameters at the first four neighbors, the second- 
order corrections to them on a ferromagnetic background, correction to the energy E_, the 

supcrcxchaugx! constant J, and the ratio li/J for the upccial case at U,/ = 8.13 o\ . 



onc-holc local 


E/ = -6.42 


Eg = 0.31 


Eg = 1.02 




energies, (eV) 










two-hole local 


E_ = -9.81 


Eo = -2.04 


E+ = 1.66 


E^ = -6.10 


energies, (eV) 


E_ = -5.40 


E+ = 1.33 


E^ = 2.04 




neighbor 


direct hopping 


contribution of 


contribution of 


second-order 


number, 


parameters. 


direct 0-0 hopping, 


Cu-0 hopping. 


corrections. 


n 


tn, (eV) 


(eV) 


t?f, (eV) 


Stn, (%) 


1 


0.53 


0.19 


0.34 


2.5 


2 


0.028 


0.085 


-0.057 


79.5 


3 


-0.076 


-0.043 


-0.033 


22.9 


4 


0.0056 


-0.011 


0.0166 


25.2 


correction to 




super exchange 




ratio 


the energy E_, (%) 




constant J, (meV) 




\ti/J\ 


4.66 




126.2 




4.18 
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TABLE V. 

Energies of the local states, the hopping parameters at the first four neighbors, the second- 
order corrections to them on a ferromagnetic background, correction to the energy E_, the 

suix^'cxchHUgc coiistaui J, and the ratio /i/J for IIk^ gcvurral case al U,/ = cS cV, A = 5.1 (A\ 



onc-holc local 


E/ = -6.43 


Eg = 0.31 


Eg = 1.02 




energies, (eV) 










two-hole local 


E_ = -9.85 


Eo = -2.12 


E+ = 1.61 


E^ = -6.12 


energies, (eV) 


= -5.41 


E+ = 1.33 


E^ = 2.04 




neighbor 


direct hopping 


contribution of 


contribution of 


second-order 


number, 


parameters. 


direct 0-0 hopping, 


Cu-0 hopping. 


corrections. 


n 


tn, (eV) 


(eV) 


t?f, (eV) 


Stn, (%) 


1 


0.528 


0.187 


0.341 


2.5 


2 


0.0274 


0.0846 


-0.0572 


82.2 


3 


-0.0758 


-0.0424 


-0.0334 


23.0 


4 


0.0057 


-0.011 


0.0167 


24.8 


correction to 




super exchange 




ratio 


the energy E_, (%) 




constant J, (meV) 




\ti/J\ 


3.02 




126.4 




4.18 
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TABLE VI. 

Contributions of p — d hopping, p — p hopping, Yp^ and Up terms to the effective hopping 
parameter of the lowest singlet at the nearest neighbors, the superexchange constant J and the 
ratio t/J at Ud = 8 eV, Up = 4 eV, Vpd = 1.2 eV, A = 3.0 eV. 



effective hopping 
parameter, 
tu (eV) 


p — d 
contribution, 
(eV) 


p — p 
contribution, 
(eV) 


contribution, 
(eV) 


Up 

contribution, 
(eV) 


0.427 


0.351 


0.202 


-0.113 


-0.013 




superexchange 
constant J, (meV) 




ratio 
\ti/J\ 






126.26 




3.38 
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